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Abstract. A formula expressing free cumulants in terms of Jacobi parameters of the 
' corresponding orthogonal polynomials is derived. It combines Flajolet's theory of con- 

tinued fractions and the Lagrange inversion formula. For the converse we discuss Gessel- 
Viennot theory to express Hankel determinants in terms of various cumulants. 
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1. Introduction 



In his seminal paper || Flajolet gave a combinatorial interpretation of continued frac- 
tions and derived generating function expressions for enumeration problems of various 
lattice paths. This is connected to work of Karlin and McGregor and others on birth and 
death processes, see [Hi] or [H| for a survey. The basic principle has been rediscovered 



many times, see e.g. |TJJ for an application to random walks and [IJ for an interpretation 
in noncommutative probability. Since Stieltjes' times continued fractions have also been 
a basic ingredient of the theory of orthogonal polynomials. A synthesis of both aspects 



. can be found in Viennot's memoir p7| or in the survey pq 



Orthogonal polynomials and their Jacobi operators have been extensively studied in 
connection with moment problems and spectral theory 0. Spectral theory of convolution 
operators was also one of Voiculescu's movitations for the development of free probabil- 

O " 



ity theory |29| and it would be interesting to understand the behaviour of orthogonal 



polynomials under free convolution. The special case of free projections, which is up to a 
translation equivalent to the case of free generators, have been considered in the literature, 
^| see e.g. H,|,|2|,|. 

In this note we derive a formula for free cumulants in the spirit Flajolet's and Viennot's 
theory of lattice paths. After discussions with A. Lascoux, he gave a proof in terms of 
symmetric functions of a more general formula for the coefficients of powers of continued 



fractions, cf. 16, 35 



The paper is organized as follows. 

In section |2] we review Flajolet's formula for the generating function of Motzkin paths. 
In section |] we briefly discuss the relevant facts about orthogonal polynomials. 
In section [| we review the definitions of various cumulants. 

In section [5] we prove the main result, a formula expressing free cumulants in terms of 
Jacobi parameters. 

Finally in section ||| we indicate how Gessel-Viennot theory of Hankel determinants can 
be used to express the Jacobi parameters in terms of cumulants. 

2. Enumeration of Lattice Paths 

Definition 2.1. A lattice path is a sequence of points in the integer lattice Z 2 . A pair 
of consecutive points is called a step of the path. A valuation is a function on the set of 
possible steps Z 2 x Z 2 . A valuation of a path is the product of the valuations of its steps. 
In the rest of this paper all lattice paths will have the following properties. 
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1. Starting point and end point lie on the x-axis. 

2. The ^-coordinates of all points are nonnegative. 

3. In each step, the x-coordinate is incremented by one. 

Thus a path of length n will start at some point (xo,0) and end at (xq + n, 0). The 
valuations will be independent of the x-coordinates of the points. Therefore the x- 
coordinates are redundant and we can represent a path tt by the sequence of its y- 
coordinates (7r(0), 7r(l), . . . , 7r(n)). We call a path irreducible if it does not touch the 
x-axis except at the start and at the end. Every path has a unique factorization into 
irreducible ones. 

We will be concerned with two types of paths. 

Definition 2.2. (i) A Motzkin path of length n is a lattice path starting at Pq = (0, 0) 
and ending at (n, 0), all of whose y coordinates are nonnegative and whose steps are 
of the three following types. 
rising step: (1,1) 
horizontal step: (1,0) 
falling step: (1,-1) 

Motzkin paths without horizontal steps are called Dyck paths. We will denote the 
set of Motzkin paths by Ai n and the subset of irreducible Motzkin paths by 
Motzkin paths are counted by the well known Motzkin numbers. 
(ii) A Lukasiewicz path of length n is a path starting at (0, 0) and ending at (n, 0) whose 
steps are of the following types. 

rising step: (1,1) 

horizontal step: (1,0) 

falling steps: (1, —k), k > 
We denote the set of Lukasiewicz paths of length n by C n and the subset of irreducible 
Lukasiewicz paths by C 1 ". Lukasiewicz paths form a Catalan family, i.e., they are 
counted by the Catalan numbers. 

The name Lukasiewicz path is motivated by the natural bijection to the Lukasiewicz 
language, the language understood by calculators using reverse polish notation, which was 



used by Raney in |23| to give a combinatorial proof of the Lagrange inversion formula. 

Flajolet's formula expresses the generating function of weighted Motzkin paths as a 
continued fraction. 

Theorem 2.3 ([g). Let 

(2.1) fj, n = 

TTeMn 

where the sum is over the set of Motzkin paths tt = (tt(0) . . . 7r(n)) of length n. Here 

is the level after the j-th step, and the valuation of a path is the product of the valuations 

of its steps v(tt) = Yli v i> ^ e latter being 

1 if the i-th step rises 

(2.2) Vi = v(rc(i — l),n(i)) = t a^j-i) if the i-th step is horizontal 



Then the generating function 



^Tr(j-i) if the i-th step falls 



n=0 
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(2.3) 



MO) 



1 — a z — 



X 1 z' 



1 — OL\Z 



X 2 Z 2 



3. Orthogonal Polynomials 

A sequence of (formal) orthogonal polynomials is a sequence of monic polynomials P n (x) 
of degree deg P n = n together with a linear functional \x on the space of polynomials with 
moments fi(x n ) = [i n such that fi{P m {x) P n (x)) = S mn s n for some coefficients s n . We will 
always assume that so = 1, that is fi(x°) = 1. It is then easy to see that such a sequence 
satisfies a three-term recurrence relation 

xP n (x) = P n+ i(x) + a n P n (x) + A„P n _i(x) 

with Ao = 0. The numbers a n , X n are the so-called Jacobi parameters. One verifies easily 
that the polynomials P n (x) = satisfy the orthogonality condition, where 



(3.1) 



and 



(3.2) 



Djx) 



Ml M2 



Mn-1 
1 X 



Mo Ml 
Mi M2 

Mn Mn+l 



M« 

M"+l 
M2n-1 



Mn 
Mn+1 

M2n 



is the n-th Hankel determinant. 

Coming from an operator theoretic background, the most natural way to make the 
connection to Flajolet's formula is perhaps via matrices. The Jacobi matrix model for 
the moment functional is 

"«o 1 

Ai «i 1 

A 2 a 2 1 



J 



that is, if we denote the basis of the vector space by e ,ei,..., with inner product 
(e m5 e n ) = S mn s n , then Je n = e n+ i + a n e n + A n e n _i. Then it is easy to see that 
P n {J) Co = e n and therefore J satisfies (J n eo,e ) = Mn- Expanding the matrix power 
yields 

Mn ^ , JQ,i\Ji\,i2 • • • Ji n -±fi 

il,...,i„_l>0 

and because of tridiagonality the sum is restricted to indices — ij\ < 1. The sum- 
mands Jo^J^,^ ■ ■ ■ Ji n _ u o can be interpreted as the valuations of the Motzkin paths 
(0, H,i 2 , ■ ■ ■ , i n -i, 0) with weights v(y, y') = J y+ iy+i, where J^-i = X { , J iti = a u J i>i+1 = 
1 and this is exactly Flajolet's formula ( p.l|) . 
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4. CUMULANTS 

Cumulants linearize convolution of probability measures coming from various notions 
of independence. 



Definition 4.1 (f29fl). A non- commutative probability space is pair (.4, ip) of a (complex) 
unital algebra A and a unital linear functional <p. The elements of A are called (non- 
commutative) random variables. The collection of moments /i n (a) = V 9 ( a ") °f such a 
random variable a & A will be called its distribution and denoted fj, a = (fi n (a)) n . 

Thus noncommutative probability follows the general "quantum" philosophy of replac- 
ing function algebras by noncommutative algebras. We will review several notions of 
independence below. Convolution is defined as follows. Let a and b be "independent" 
random variables. Then the convolution of the distributions of a and b is defined to be 
the distribution of the sum a + b. In all the examples below, the distribution of the sum 
of "independent" random variables only depends on the individual distributions of the 
summands and therefore convolution is well defined and the n-th moment /i n (a + b) is 
a polynomial function of the moments of a and b of order less or equal to n. For our 
purposes it is sufficient to axiomatize cumulants as follows. 

Definition 4.2. Given a notion of independence on a noncommutative probability space 
(A,<p), a sequence of maps a t— > k n (a), n = 1,2, . . . is called a cumulant sequence if it 
satisfies 

1. additivity: if a and b are "independent" random variables, then k n (a + b) = k n (a) + 
k n (b). 

2. homogenity: k n (Xa) = \ n k n (a). 

3. k n (a) is a polynomial in the first n moments of a with leading term /i n (a). This 
ensures that conversely the moments can be recovered from the cumulants. 

We will review here free, classical and boolean cumulants via their matrix models. The 



reader interested in "g-analogues" is referred to the g-Toeplitz matrix models of jl9] and 
which yield similar valuations on Lukasiewicz paths. 



4.1. Free Cumulants. Free probability was introduced by Voiculescu in |ID| and has 
seen rapid development since, see and the more recent survey 



Definition 4.3. Given a noncommutative probability space (A, <p), the subalgebras A% C 
A are called free independent (or free for short) if 

(4.1) p{aia 2 ■ ■ -a n ) = 

whenever aj G Ai- with ip(a,j) = and ij ^ ij + i for j = 1, . . . , n — 1. Elements dj G A are 
called free if the unital subalgebras generated by a, are free. 

Existence of free cumulants with the properties of definition |4.2| was proved already 
in [30]. A beautiful systematic theory of free cumulants was found by Speicher in his 



combinatorial approach to free probability via non-crossing partitions |24 |; see [|17]] for 



an explanation why noncrossing partitions appear. The explicit computation involves 
generating functions as follows. 

Theorem 4.4 (0]). Let M(z) = 1 + Y^Li ^n- 2 ™ be the ordinary moment generating 
function and C(z) = 1 + Y^=i c n zn be the function implicitly defined by the relation 
C(zM(z)) = M(z). Then the coefficients c n satisfy the requirements of definition 
and are called the free or non-crossing cumulants and C(z) is the cumulant generating 
function. 

The moments and the free cumulants are related explicitly by the following combinato- 



rial formula (see [^3|, § for the formulation in terms of Lukasiewicz language, and |24| for 
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non-crossing partitions). Define a valuation on Lukasiewicz paths by putting the following 
weights on the steps 

v(y, y — k) = c k+ i k > 

(4.2) ) + 

then 

(4.3) Mn =5>(7r). 

Another interpretation of this relation is the Fock space model of Voiculescu (see 
13| for a simpler proof). It involves Toeplitz operators as follows. Let S be the forward 
shift on ^2 (No), i.e., Se n = e n+ \. A (formal) Toeplitz operator is a linear combination 
of powers of S and of its adjoint, the backward shift S*e n = e n _i. The linear functional 
u)(X) = (Xeo,eo) is called the vacuum expectation. Let c n be the free cumulants of the 
moment sequence (fi n ) as defined above and set 

oo 

(4.4) T = S* + Y,c n+ iS n . 

n=0 

Then u(T n ) = fi n and writing T in matrix form 

"ci 1 

T 



c 2 Ci 
c 3 c 2 



1 



we can expand the matrix product and obtain 
u(T n ) = (T n e ,e } = 



.,i n _i>0 



To i-, Tj 



l- t *l,«2 



T, 



in-1,0 



where Tij ^ only for j > i — 1. Again this can be interpreted as a sum over lattice 



paths which this time turn out to be Lukasiewicz paths with weights (|4.2j ) coming from 
the matrix entries 



and we obtain the sum ( |4.3| ). There are multivariate generalizations of this Toeplitz 
model, see [2T | . 

4.2. Classical Cumulants. Classical cumulants linearize convolution of measures and 
can be defined via the Fourier transform. 



x Hn z n De the exponential moment generating 



Definition 4.5. Let F(z) = fi(e xz ) = J2 n=0 
function. Let K(z) = log F(z) = Yl^Li ^f 2 ™ be its formal logarithm. The coefficients n n 
are the classical cumulants of the functional \x. 

There is a model on bosonic Fock space for classical cumulants which is analogous to 
the Toeplitz model ( |4.4| ). 

Let D, x be annihilation and creation operators which satisfy the canonical commutation 
relations (CCR): [D,x] = 1, i.e., on the Hilbert space with basis e n and inner product 
(e m , e n ) = n\5 m>n , set De n = ne n _i, xe n = e n+ i. Then D and x are adjoints of each other. 

Denote again by u> the vacuum expectation u(T) = (Teo, eo), i.e. u(x k D n ) = Wn ^ 
and u(f(x)) = /(0). Then the Fourier-Laplace transform of 

oo 

l~i I, 



T = D + y^+i x n 



is 



io[e 



n=0 



, n — 1 
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In the basis {e n }, T has matrix representation 



T 



'k\ 1 

K 2 K\ 2 

^ K 2 Ki 3 

HA HI fc„ AC, 

31 2! 2 1 



e.g. for the stnadard Gaussian distribution we have n 2 
so the model is tridiagonal 

"0 1 



1 and higher cumulants vanish, 



1 
1 



2 




and this is the well known Jacobi operator matrix for the Hermite polynomials. Similar 
to the formula above we get the following combinatorial sum for the moments of T: 



with valuation 
(4.5) 



v(y,y- k) 



k > 



v(y,y+ l)=y+l 

4.3. Boolean Cumulants. Boolean cumulants linearize Boolean convolution which was 



introduced in [^IJ (compare also [33], [[34]) and is a special case of conditional indepen- 
dence (0, @])- In the theory of random walks boolean cumulants arise as first return 
probabilities. 

Definition 4.6. Given a noncommutative probability space (.4, ip), the subalgebras Ai C 
A are called boolean independent if 

(4.6) ¥?(aia 2 • • • a n ) = yj(ai) ^(a 2 ) • • ■ y?(a n ) 

whenever a,- G 4^ with 7^ ij + \ for j = 1, . . . , n — 1. 

With this notion of independence, convolution of measures is well defined and the 
appropriate cumulants can be calculated as follows P5[] . Let M(z) = X^o/- 4 ™- 2 ™ ^ e ^ ne 
ordinary moment generating function. Then 



(4.7) 

where 

(4.8) 



M(z) 



H(z) 



1 - H(z) 



n=l 



where h n are the boolean cumulants of the distribution with moment sequence (/i n ) . Ex- 
panding (|4.7|) in a geometric series we find that the moments can be expressed as Cauchy 
convolution 



E E 

r=l ijH h*r=™ 

i,->l 
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and from Theorem |2]3| and the fact that every Motzkin path has a unique factorization 
into irreducible paths, it follows that 



(4.9) h n = 



where v is the valuation Q2.2|) . This can also be seen by comparing the continued fraction 
expansion ( [2.3|) with (|4.7| ). In terms of the Jacobi parameters the boolean cumulant 
generating function has the continued fraction expansion 

H(z) = a z H — 2 ■ 



l- aiZ 

1 — a 2 z 



Similarly, we can express the boolean cumulants in terms of free cumulants (resp. 
classical cumulants) 



(4.10) h n = 
using the valuations (|4.2|) (resp. (|4.5|)). 



5. A Formula for Free Cumulants 

Theorem 5.1. For n > 2, we have the following formula for the free cumulant c n in 
terms of the Jacobi parameters: 

(kk yr (-i) Ho ~ x fn-i\ 

(5 - 1} Cn = 22 -i^rr{\«\ ) v W- 

TteMn \ I 10 / 

Here \tt\ is the number of returns to zero of the path and v is the valuation fl2.2p from 
Flajolet's formula. Note that the path consisting of n horizontal steps does not contribute. 

Proof. The proof consists in comparing the expressions for the moments and free cumu- 
lants in terms of the Boolean cumulants (|4.8|) , which themselves represent the sum over 
irreducible Motzkin paths ( |4.9| ). We have then by the multinomial formula 

Hn = [z n }M(z) 
n 

= [z n ] H{z) m 



m=0 

n 

+2fc2H \-nk n • 

ii/. ''ii 

k\ , . . . , k n 



m=0 fci+fe2H \-k n =m 



E (\ + '"t fc, )*f-^ 



. k\ ..... kn 

ki~\-2k2~\ \-nk n =n N 
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The cumulants can be expressed in terms of the moments with the help of Lagrange's 
inversion formula 0: 

L J 71 JT l -An.— 1 



1 


n 




1 




1 




n 




i 
i 




1 




n 




1 




1 




n 




1 




1 




n 




1 



z n }(l- H(z)) n - 1 

z n ](l -h x z Kz n ) n ~ l 

z n \ ^~ ^ ! , (-/ii) fcl • ■ • (-h n ) kn z^ +2k2+ - +nkn 



k \ ■ ■ ■ k n \ 

fco+A:iH hfc„=n-l 



fei+2A;2H \-nk n =n 

ki <n 



— 1 . , \n 1 k\ ••• /c n , k± } . . . j fi n / 



^ n-l Ui + --- + W\ h,...,k n J 1 



n-l , , w„! 



e e ^r;V-^ 



n — 1 \ r 

r=l iiH Yi r =n x 



Comparison with the formula for the moments yields the result. □ 



Remark 5.2. Actually by the same argument we can express free cumulants in terms of 
classical cumulants by simply replacing Motzkin paths by Lukasiewicz paths with valua- 
tion ( |4.5|) . However, there are cancellations in the sum. There are also lots of cancellations 
when expressing free cumulants in terms of free cumulants themselves using Lukasiewicz 
paths with valuation Q4.2| ) (only one term survives). Note however that there are no can- 
cellations in the sum Q5.1J) , since \tt\ is equal to the number of ao's and Ai's and therefore 
every monomial appears always with the same sign. 



6. Hankel determinants 



In this section we survey some formulas expressing the Jacobi parameters a n and \ n in 
terms of cumulants. 

It is well known (see e.g. @) and can be readily deduced from ( |3.1|) , that the Jacobi 
parameters can be expressed in terms of Hankel determinants and Hankel minors of the 
moments as follows: 

A n _2A n 



A, 



A 2 



and 



where 



A, 
A r 



A 



n— 1 



A 



A r 



//•o 
/'i 



/<2 



n— 1 



fJ-n fJ'n+l ■ ■ ■ ^2n-l /^2n+l 

We are therefore interested to express the minors of the infinite Hankel matrix H,^ 



i+j\i,j>0 



in terms of cumulants. Namely, given finite sequences of indices ii,i2, 
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and ji, j 2 , . . . , jp we will denote 



(6.1) 

For example, 

A„ = 



H 



n 12 



<"p 
h 



A^2 + jp 



H 



1 
1 



A n = H 



1 
1 



n-1 



n 



n 



1 n + 1 



The main tool is the so-called Gessel-Viennot theory, see |27| [12|, [15], [18|]. Let T be 
a weighted graph, that is a graph with vertices V and edges E together with a val- 
uation v : E — > on the edges, where Z? is some (commutative) ring. The valua- 
tion of a path cj = (s , • • • , s m ) is the product of the valuations of the steps v(u) = 
v(s , si) v(s u s 2 ) ■■■ u(s m _i, s m ). 

Select two sets of distinct vertices Ai, Bi G V, i = 1, . . . , n and suppose that the set of 
paths 



{u = (s Q , . . . , s m ) : s = A, s m = 5,-, v(u) ^ 0} 



is finite for every pair of indices We will say that two paths intersect if they have a 

vertex in common. We will not care about crossings of paths, that is, crossings of edges 
in a graphical representation. 
Define a matrix 

Then the determinant of this matrix has the following combinatorial interpretation. 
Proposition 6.1 (071 Prop. IV.2). 



(6.2) 



det[a i:; -] nxn = 22 siga(cr) v(ux) • • • v(u n ) 

(a;u>l,...,u>„) 



where the sum is over all permutations a G & n , iOi G f2j )CTj and all non-intersecting (but 
possibly crossing) paths u>i G £\o-(i)- 

We have seen various examples above where the n-th moment is equal to the sum of 
the valuations of paths of length n. All our valuations are translation independent and we 
will interpret the entry of 7i M as a sum = J2 V ( U} ) over the set of paths starting 
at (— i, 0) and ending at (j, 0). The minor (|6.1| ) is therefore equal to 



H 



ll 12 

3i k 



Jp 



E 



E < 



where the sum extends over all permutations a of the indices and all nonintersecting 
configurations of paths Uk starting at (—ik, 0) and ending at (j a (k), 0). Because the vertices 
are ordered, the sign is (— l) sl s n ( 7r ) = (— 1) K , where K is the number of crossings of the 
configuration corresponding to n. 

For example, for the full Hankel determinant from (|3.2|) , where the moments are inter- 
preted as sum over Motzkin paths as in Theorem |2.3| , there is only one non-intersecting 
configuration whose picture is (for n — 4) 
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H 



and therefore A n = A™A^ 1 • • • X^_ x \ n . 



0,1,2,3,4^ 
0,1,2,3,4. 

A n- 

Note that in this example the point (0,0) is considered as a path of length zero and is 
not allowed to lie on any of the other paths. 

In contrast, in the case of Lukasiewicz paths with valuations (|4.2| ) and (|4.5| ), there may 
occur crossings of paths and there are many contributing terms with different signs. It is 
also possible that cancellations occur in the sum, e.g. 

/ \ / \ 




V 



are cancelling contributions to H 



2 2 

-c x c 2 c z 




\ 



J 



1 2 
1 3 



. Therefore, and because the sums are com- 



plicated, the formulae are of rather limited value. 
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